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Abstract 

In this paper we have explored the consequences of a model of 
dark energy with its energy density varying exponentially with the 
scale factor. We first consider the model with p^ oc e Ka , where k is a 
constant. This is a kind of generalisation of the cosmological constant 
model with k = 0. We show that such an exponentially varying dark 
energy density with the scale factor naturally leads to an equivalent 
phantom field. We also consider a model with p^ oc e K / a and we show 
that this also naturally leads to an equivalent phantom field. 
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1 Introduction 



Observations on Supernova Type la, CMB and studies on large scale struc- 
ture formation have strongly suggested that at the present epoch the universe 
is in a phase of accelerated expansion P3I2]- Within the standard framework 
of General Relativity, this kind of expansion is not possible with the con- 
ventional kind of source, i.e., which has an equation of state of the form 
p = wp with < w < 1. The present accelerated expansion of the universe 
seems to indicate that there exists some form of matter-energy with negative 
pressure which is causing this acceleration. (One could, of course, argue that 
Einstein's general relativity may not be the fully correct /complete theory on 
the cosmological scales |3 . But we will not consider this possibility in this 
paper.) The former possibility has motivated the exploration of alternative 
kinds of sources with p < — |p as well as a resurgence of investigations us- 
ing cosmological constant (w = —1). Such a matter energy component with 
equation of state p < — |p is generally known as 'Dark energy'. Several dark 
energy models have been considered in literature. The cosmological constant 
is just the simplest case with w — — 1. Models with cosmological constant 
predict a transition to an accelerated expansion phase of the Universe. 

Introducing the cosmological constant, however, leads to two serious prob- 
lems. The first one is that the observations indicate that the value of A is 
many orders of magnitude smaller than what is expected from the energy 
scales involved in the Early Universe. This problem is known as the cosmo- 
logical constant problem [3J |H| • The second problem is that the observa- 
tions indicate that at the present epoch the energy density of the A term, 
i.e., pa — is of the same order in magnitude as the energy density of 
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matter p m with A having begun to dominate in the recent past. This re- 
quires extreme fine tuning of the value of A. This problem is known as the 
cosmic coincidence problem. The vacuum energy density p\ does not change 
with the expansion of the Universe . Hence, there is a need to explore more 
general models of dark energy. 

The other models of dark energy in which the density of dark energy, pd e , 
changes as the universe expands are based on Quintessence, Phantom field, 
Tachyon field etc. Quintessence models refer to those models in which dark 
energy is represented by a standard scalar field with a potential V(4>) [HEHH]- 
Specific forms of potential which satisfy certain conditions admit tracking 
solution which may possibly solve the cosmic coincidence problem ^UJ ■ 
Quintessence models have an equation of state parameter u^ e greater than 
minus one. Observational constraints on Wd e do not exclude the possibility 
of Wde being less than minus one at the present epoch jT2l EE] • This leads 
to the possibility that the present cosmic acceleration may be driven by a 
phantom field (wd e < — 1). Energy density of the phantom field increases as 
the universe expands. The phantom energy could be obtained from either 
a scalar field with negative kinetic energy ^] term or by some form of 
non minimal coupling of the scalar field with gravity ^B]. The phantom field 
described in this paper refer to that with a negative kinetic energy term. 
Some models of phantom energy with a constant Wd e lead to a singularity 
after a finite time in the future [T71 118j. This is called the Big Rip. Certain 
other models with dynamical Wd e do not lead to Big Rip There is yet 
another theoretical possibility that the present accelerated expansion might 
be due to tachyon fields EI] ■ These were originally motivated from the 
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string theory. In these models Wd e > — 1- 

In all these scalar field models the equation of state is given as p$ = w^p^,. 
Dark energy with different equation of state has also been considered in 
literature, for example, Chaplygin gas with p = - [211 121] and models of 
dark energy with generalised equation of state p = a(p — p Q ) [22] . 

It has been argued that probing the dark energy via its density evolution 
rather than its equation of state has a number of advantages [213 • I n this 
paper we consider a model of dark energy with its energy density varying 
exponentially with the scale factor. After discussing the basic formalism for 
a general scalar field in section [21 we will consider in section El the model 
with oc e Ka where k is a constant. This is a kind of generalisation of the 
cosmological constant model which corresponds to k — 0. We will show that 
such a model naturally leads to an equivalent phantom field. This model, 
however, leads to Big Rip. This motivated us to consider another model with 

oc e K / a which is discussed in the section 0] This also naturally leads to an 
equivalent phantom field but it doesn't lead to Big Rip. It asymptotically 
behaves as a de Sitter Universe. The phase space analysis has been discussed 
in sectional 

2 Energy- momentum Tensor and the Equa- 
tions of Motion for a general scalar field 

Let the Lagrangian density of the scalar field be given by and that for 
the rest of the matter by, C source . The complete Einstein-Hilbert action for 
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this system is then given by, 

-1 



S 



lo7rG 



"=9 d A x (1) 



We employ the metric with signature (H ). The Lagrangian density 

for a normal scalar field is given by, 

c^h^d^-v^) (2) 

The kinetic energy term for this scalar field is positive. Another class of 
scalar fields have been considered in literature for which the kinetic energy 
term is a negative. These are called phantom fields and their Lagrangian 
density is given by, 

C^-U^d^-V^) (3) 

In order to incorporate both these kinds of scalar fields, we consider an action 
of the form, 

S = J {^Q R + [f ^ W " W)] + £source}V=gd 4 X (4) 

where a takes the value +1 for normal scalar field (quintessence) and —1 
for phantom field. The energy momentum tensor of the scalar field obtained 
from the action in equation (@J is given by: 

2$ = ad^d^ - gr^d,^* ~ V(<P)] (5) 

Consider the spatially flat FRW metric: 

ds 2 = dt 2 - a(tf[dr 2 + r 2 d6 2 + r 2 sin 6 2 d<f) 2 ] (6) 

We further assume that the scalar field is homogeneous. The spatial 
derivatives of such a field vanishes and the field depends only on time i.e 



= 4>(t). Using this form of and the metric corresponding to the line 
element in equation ©, the energy momentum tensor in equation (jSJ) assumes 
a diagonal form given by: 

= diag(p <i> , -P4>, -P<j» ~P<t>), (7) 

where, 

P<t> = «y + V{4>) (8) 



and 



«y - (9) 



The equation of state is given by = w^p^ where, 

gj[ - ^(0) nm 

^ = — n ( 10 ) 

«f + U(0) 

The dynamics of the Universe is determined by the Friedmann equations: 

H* = 8 -fp (11) 
^ = -^(P + 3 P ) (12) 

where p and p are the total energy density and the total pressure. ( We have 
assumed a k = universe). The field equation for is obtained from the 
conservation equation T% v . v = and is given by: 

'(p + m^ + a' 1 ^- = (13) 
dip 

Let p = p m + P<t> (We have neglected the contribution from radiation since 
we are interested only in the late time evolution of the Universe). The scale 
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factor is normalised so as to make a (to) — 1- The density of matter is hence 
given by 

Pm = Pmo<^ ?J (14) 

where p mo is the matter density at the present epoch. 



3 Model with oc e Ka 



As we stated in the introduction, our approach in this paper is to start with 
a form for the evolution of the energy density of the scalar field rather than 
the equation of state. This is the approach that has been adopted in [25] . 
The present paper is in this spirit. We consider a scalar field in a potential 
V(4>) whose energy density varies with scale factor as p^ oc e Ka . Since we 
have normalised the scale factor to unity today, we may write, 

P<t> = P0oe K{a - 1} (15) 

where p^o is the dark energy density in the present epoch with the scale 
factor today normalised to unity. The continuity equation for the scalar field 
T£ v . v = implies: 

3 

— (l + w^pcj, (16) 



da a 

where, is the ratio of the pressure, p^, to the energy density, p^, of the 
scalar field. Using the form for as given in equation (jl5j) . we get, 

w+ = -\-— (17) 

From equation (j!7|) it follows that if k < then at sufficiently large value 
of the scale factor a , the ratio of the pressure to density of the scalar field, 



would become greater than one. This would violate causality. In order 
to ensure causality we need < +1 and in order to ensure this at all future 
epochs, we require k to be greater than zero. 

Using equations (JBJ) and (JTUJ), we have 

a<j) 2 = (l + w^ P(j> (18) 
Substituting for from equation (fT7j) . we get, 

= - yP* (19) 

For > and k > (as required by the causality condition ), a should be 
negative. Hence, the kinetic energy term must have a negative sign which 
fixes a = — 1. Hence such a form of p^ naturally leads to an equivalent 
phantom scalar field. 

We now explore the nature of the potential V((f>) of this phantom field 
which would lead to such a form of p^. From equations (jHJ), (jUJ) and (fTUj) 
(with a = — 1) we have, 

P<p = -y+V(0 (20) 
= ~ -V(0) (21) 

= 4^ (22) 



-f + ^(0) 



(23) 



These equations imply, 



2 = -(1 + W ^ (24) 
7(0) = (25) 



The above equation (|25jl together with equations (|TS|) and (fT7j) gives the 
potential V(<p) as a function of the scale factor a. The derivative of (f> with 
respect to scale factor a is given by: 

f = t = (26) 
da d aif 

where if is given by the Friedmann equation (fTTj) : 

^^(p^ + P^ 1 ^) (27) 

Using equations JHJ), (jnj), (J23), (ESJ), flU and flZZJ) one obtains, 
d <j) a(zd^) 1 / 2 



(28) 



Here, 



and 



And 



V(a) = ^HlMlil + ^e^. (29) 



n rn = p r - (31) 



M 2 = — (32) 



which is the reduced Planck mass in units of % = c = 1 . -£/o is the Hubble 
parameter. Using equations (|28|l & ()29|) one can plot potential V(<p) as a 
function of 0. This plot is shown in figure (|T|). 

In this plot we have chosen initial condition such that <p = at a = 10 -5 
. The field (f> grows from (f) = to a value = 0.55M p at the present epoch 
(a=l). 
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Figure 1: This figure shows the potential V((f)) for the model oc e . In 
this figure X = and Y = ^-ifji 10~ 4 . For this plot we have taken k — 1, 

pi pi 



= 0.73 and fL 



0.27. 



At low value of the field i.e. for <p M p which occurs at a < 1, the 



differential equation (J2SJ) can be approximated as 



d (f) I kVL 

— ( ) ~ a 

da M„ 



^ Lr m x ~' 



(33) 



which on integration gives: 



/ a 



(34) 



Where O is the constant of integration and 



a = 2 A 



/ O m 6 



(35) 



On substituting equation (J3"4*j) in the equation (J2llj) . we obtain an approximate 
functional form of V(<j)) for small value of the field (<j) -C M p ). This is given 
as: 



Vtf) = V (l + - <Po))exp(J ^ - O )), 



(36) 
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where 

Vq = 3fi^ 2 M p 2 e - K (37) 

Similarly for large value of the field, i.e. for <p > M p , which occurs at 
a > 1 (and hence at a future epoch), the differential equation (J25J) can be 
approximated as: 

which on integration gives: 

- ( ^« )2 (39) 



AkM% 

Here O is the constant of integration. Substituting equation in the 
equation (|29|). we obtain an approximate functional form of the potential 
V{4>) at large value of the field <fi. This is given as: 

V{4>) = V (l + - 10? -) e x P (_ 0? _) (40) 

We will consider this form of potential given by equation ()4Uj) in section El to 
study the phase space analysis. 

The value of the constant k can be determined from the epoch of equality 
of matter and dark energy density. Let z m ^ be the red shift of the epoch 
of matter-dark energy equality. Equating the expression for and p m from 
equations (fT5|) and (jTlj) we get, 

« = -(l + — )ln(^(l + z m(f> f) (41) 



As k > equation (PT|) implies, 



W < (^) V3 - 1 (42) 
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As = 0.73 and Vt m = 0.27 (as indicated by the observations by WMAP 
|26j) we have z m ^ < 0.393. This model (with p^ oc e Ka ), thus, constrains the 
value of We can see from the expression of w<j>(a) given by equation lfTTj) 
that a — > oo =>- w</,(a) — > — oo. From the expression of p^ (fT5j) and the 
equation of state we see that p^ — > oo and — > — oo as a — > oo. This would 
then mean that the trace of the energy momentum tensor T = Tjf = p^ — 3^^ 
(and hence the Ricci scalar R oc T) would blow up as a — >• oo. It is of interest 
to ascertain the epoch when the scale factor tends to infinity. In particular, 
we would like to know whether this happens at a finite time or asymptotically 
in the infinite future. 

The epoch, t^, when the scale factor becomes infinitely large can be 
calculated from the Friedman equation (|27jl to be, 



large, the first term under the square root term becomes subdominant as 
compared to the second and so the integral remains finite. Thus this model 
exhibits a singularity at a finite future epoch. This is popularly known as 
Big Rip and our model p^ oc e Ka exhibits that feature. 




(43) 



where too is the epoch when a 



oo. As the scale factor becomes infinitely 



4 Model with oc e 



k/cl 



Let p = Pm + P<t> where p^ oc e K / a . Denoting by p^o the density of the scalar 
field in the present epoch corresponding to a = 1, we have: 



X-a 
K 



(44) 
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We again consider the action of the form (J1J. The continuity equation 
T^ u .„ = would imply that: 

W * = ~ 1 + Ta (45) 

If k > then at sufficiently low value of the scale factor a, would become 
greater than one. This would violate causality. So for to be less than 
one requires k to be less than zero. Hence in this model causality requires 
that k < 0. With these constraints, we investigate the sign of the kinetic 
energy term of the scalar field that would lead to p^ of the form (J44|) and 
also calculate the form of the potential V(4>) of the field 0. Using equations 
©, ©, (EH), il & (US) we obtain, 

«0 2 = ^ (46) 

For > and k < (we have shown that for this model causality requires 
k < 0), the left hand side of the equation is negative. So in this case too, the 
kinetic energy term must be negative. Hence such a form of p^ also naturally 
leads to an equivalent phantom scalar field. We now obtain the form of the 
potential V{4>) of the phantom field which would lead to such a form of p$. 
For p$ given by equation (j33j), the Friedmann equation (|27j) becomes: 

H 2 = ^-(Pmoa~ 3 + P^e^) (47) 
Using equations (JSj), (JUJ), 421, 425), 42ED and (JUD we obtain, 



(48) 
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this figure X = and Y = H -z?k ■ For this plot we have taken k = — 1, 

pi pi 

fi = 0.73 and fi m = 0.27. 

Using equations (JlT^) & (J4T?|) one can plot potential V(4>) as a function of 0. 
This plot is shown in figure @. In this plot also we have chosen an initial 
condition such that = at a = 10~ 5 . This gives a value of the field 
<p = 0.7M P at the present epoch (a = 1). The field grows from </> = to 
a value <p ~ 2.45M P . As done in the previous section we obtain the limiting 
form of the potential V(4>) from equations (}4*Sj) and (|4*U|) . 

When <C M p which occurs at a -C 1, the potential V(4>) is nearly flat. 
When cf) > M p and a ^> 1, the equation can be approximated as: 



which on integration gives: 



Here 0o is the constant of integration. Substituting this equation in equation 
(149)1 . we obtain a limiting form of the potential at large value of the scale 




(50) 
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factor and for > M p . This is given as: 



where 

V Q = 3^iJ 2 M p 2 e- K (53) 

It is interesting to note that this limiting form of the potential given by equa- 
tion ()52)l is very similar to that obtained in the previous section ( equation 
f!40|) ) with the only difference in the sign in the exponential part. We will 
also consider this form of the potential given by equation (p)2|) in section |S] 
to study the phase space analysis. 

In this model also we can obtain the value of the constant k from the 
red shift z m ^ of the epoch of matter-dark energy equality. Equating the 
expressions for p^ and p m we obtain: 

n = —ln{^{\ + z m4> f) (54) 

As k < 0, the above equation (JSH) would imply that: 

< (§^) V3 - 1 (55) 

This is same as equation (J42)) . Hence we can see that both this model and the 
model considered in the previous section puts same constraints on z m $. Using 
equations (jBj) & (|4"5j) we can see that this model (p^ oc e K ^ a ) asymptotically 
behaves as a De Sitter's Universe. So this model does not lead to Big Rip. 
This is unlike the previous model which exhibits a Big Rip. 
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5 Phase space analysis 



As the Universe expands, the relative contribution of pressure-less matter 
keeps decreasing and as the scale factor becomes very large, the matter con- 
tent can be neglected. Let us consider the model with oc e Ka and analyse 
the situation in this limit of large a. In this limit, we neglect the energy 
density contribution of matter and only consider that due to scalar field and 
integrate the equation of motion for the scalar field to obtain. 



0(a) = 2 V8^ + ° (56) 
where 0o is the constant of integration. Using this equation and the equation 
(|29|) we obtain the following equation: 

where 

V = 3H%M*e- K (58) 

(Equation (|57j) is the same as the equation pUj) but in the equation (J57)l we 
have taken O to be zero ). From equation (J 17)) the redshift dependence of 
w^a) is of the form, = — 1 + ma where m is a constant and a is the scale 
factor, together with the continuity equation, 

^-(pa 3 ) + 3pa 2 = (59) 

this form of w^(a) imposes, p$ = Ae~ 3ma . So in order to investigate whether 
different initial condition leads to p^ of the form p^ oc e Ka , we plot as a 
function of the scale factor for different initial conditions. This plot is shown 
in fig®. From this we can conclude that for different initial condition w(a) is 
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Figure 3: Plot of equation of state parameter w(a) verses scale factor for 
different initial conditions for the potential given by equation ([57)1 . Here 
p m = and k = 2.5. 

of the form w(a) = — 1+ma asymptotically with the constant m depending on 
the initial condition. Hence, asymptotically p^ will be of the form p^ = Ae Ka 
for different initial conditions but both the constants n and A depends on 
the initial condition. Phase portrait (for the case when the potential is given 
by equation (jSTJ)) is shown in fig(HJ). 

Similarly for the model p,* oc e K / a , in the limit of large a, we get, 



V{<f>) = V (l + 



24M p 2 



)exp(- 



4M p 2 



(60) 



where 



V = ZHlM 2 p e- K (61) 



(Equation (|50jl is the same as the equation but in the equation (|50|) we 
have taken 0o to be zero ). Here we plot w(a) verses 1/a for different initial 
conditions. This plot is shown in the fig©- Here also we can conclude that 
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Figure 4: Phase portrait in model with potential given by equation (|57|). 
Here p m = and n = 2.5. The initial condition chosen is the same as in 
fig©. Here X = <p/M p and Y = <p/H M p 




-1.8 - 

-2 - \ - 

-2.2 - V 

-2.4 I 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

1/a 

Figure 5: Plot of equation of state parameter w(a) verses inverse of scale 
factor for different initial conditions for the potential given by equation (|fi0)). 
Here p m = and k = —2.5. Here 1/a — > assymptotically into the 
future. 
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-2.5 - 



-3 I 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 

X 

Figure 6: Phase portrait in model with potential given by equation (jfiOJ) . 
Here p m = and « = —2.5. The initial condition chosen is the same as in 
fig®. Here X = <p/M p and Y = <p/H M p 

for different initial conditions w(a) is of the form w(a) = — 1 + — asymptot- 
ically with the constant m depending on the initial condition. We conclude 
that asymptotically p^ will be of the form p^ = Ae K ^ a for different initial 
conditions but both the constants k and A depends on the initial condition. 
Phase portrait (for the case when the potential is given by equation (|60|)) is 
shown in fig©- 

6 Conclusions 

In this paper we have shown that the dark energy density varying exponen- 
tially as oc e Ka and p^ oc e K//a naturally lead to an equivalent phantom field. 
We have shown in section that in the limit of large a phantom potential 
of the form given by equation (|57|) leads to p^ oc e Ka for a particular initial 
condition on and 0. For other nearby initial conditions p^ tends to this 
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form asymptotically and the constant k depends on the initial condition {(f) 
and <fi). This model leads to Big Rip. Similarly in the large a limit, phantom 
potential of the form given by equation (JoUj) leads to p^ oc e K ^ a for certain 
initial conditions. For other nearby initial conditions p^ approaches the form 
oc e K / a asymptotically and the constant k depends on the initial conditions 
(0 and 0). It does not leads to Big Rip but rather it asymptotically leads to 
De Sitter Universe. 
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